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Abstract. GaN nanowire ensembles with axial In^Gai-zN multi-quantum wells 
(MQWs) were grown by molecular beam epitaxy. In a series of samples, we 
varied the In content in the MQWs from almost zero to about 20%. Within the 
nanowire ensemble, the MQWs fluctuate strongly in composition and size. Statistical 
r*j 1 information about the composition was obtained from x-ray diffraction and Raman 

O ■ spectroscopy. Photoluminescence at room temperature was obtained in the range from 

2.2 eV to 2.5 eV depending on In content. Contrary to planar MQWs, the intensity 
increases with increasing In content. We compare the observed emission energies with 
transition energies obtained from a one-dimensional model, and conclude that several 
mechanisms for carrier localization affect the luminescence of these three-dimensional 
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1. Introduction 

GaN nanowires (NWs) are pursued as an alternative to planar GaN for the fabrication 
of light emitting diodes (LEDs) [H O H U M- The advantage of GaN NWs 
is that they can be synthesized on cost-effective Si substrates in excellent crystal 
quality [7j. NWs also have the benefit over planar layers that strain from lattice- 
mismatched heterostructures can elastically relax at the free sidewalls and defects 
from plastic relaxation can be avoided [8]. Several groups have fabricated GaN NWs, 
with repeated axial Ina-Gai-^N segments to form multi-quantum- wells (MQWs), on 
Si substrates by plasma-assisted molecular beam expitaxy (PA-MBE). Luminescence 
with peak emission wavelengths from 440 nm to 640 nm has been reported [Tj, [5j [6J [9] . 
For the application of such MQWs in devices like LEDs it is mandatory to control 
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the emission wavelength. Still, at present, it is even unclear how the emission of 
ensembles of axial In^-Gai.^N/GaN NW heterostructures varies with the In content. 
The luminescence wavelength is known to increase at lower growth temperatures of the 
active NW region, because more In is incorporated [6]. However, there are only few 
reports on the determination of the In content in NW MQWs which do not rely on the 
luminescence wavelength and therefore assume its dependence on the In content a priori 
0, |T0] . Moreover, in these reports the In content has been quantified only for single 
NWs, and the effect of a systematic variation of the In content has not been shown. 
There are two major experimental difficulties. First, for NWs it is significantly more 
difficult to determine the MQW structural parameters than for planar layers. Second, 
these parameters fluctuate significantly on a given sample. Analytical tools are required 
which probe large ensembles of NWs and deliver statistical information on the MQW 
structure. Recently, we have shown that the composition and thickness fluctuations 
of In x Ga 1 _. !; N/GaN MQWs in self-induced NWs, which form superlattices if they are 
equally spaced, can be determined by laboratory x-ray diffraction (XRD) (TTJ [121 03] • 
Here, we make use of this method and analyze a series of In^Gai-^N/GaN NW samples 
for which the growth temperature of the active region and consequently the average In 
content was varied. The results for the In concentration are corroborated by resonant 
Raman spectroscopy. We then compare the observed photoluminescence (PL) with the 
expected transition energies calculated by a one-dimensional model. 

2. Sample growth 

GaN NWs were grown on Si(lll) by PA-MBE at 770 °C with a N/Ga flux ratio of 4. 
After one hour of growth, the metal flux was interrupted and the substrate temperature 
was lowered for the growth of the active region, which is formed by six identical pairs 
of GaN barriers and In^Gai-^N QWs. The achieved nominal thicknesses were 14 nm 
for the barriers and 5 nm for the QWs. During the initial phase of the growth of the 
active region, evidence for the formation of cubic GaN (c-GaN) was obtained in situ by 
reflection high-energy electron diffraction (RHEED). Figure [1] (a) shows the RHEED 
pattern during the growth of the GaN NW base, when only the reflections of hexagonal 
GaN (h-GaN) were observed. The appearance of additional c-GaN reflections after 
lowering the growth temperature is depicted in figured] (b). The inclusions of c-GaN 
affect the optical properties of the NWs, as we will see later. 

A series of five samples was produced with the growth temperature of the active 
region, T a , ranging from 583 °C to 637 °C while all fluxes and shutter opening sequences 
were kept unchanged. At these temperatures, the desorption of In cannot be neglected. 
Therefore, the resulting composition x depends on T a p3j. Scanning electron microscopy 
has revealed that the resulting NW diameters fluctuate, with typical values between 
30 nm and 100 nm (images not shown). 
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Figure 1. Evolution of RHEED pattern during growth of GaN NWs (a) in the base 
and (b) in the active region at reduced growth temperature, before In ;c Ga 1 _ a ,N is 
grown. The appearance of the cubic GaN diffraction pattern is depicted above for 
clarity (one of the twin patterns drawn). 
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Figure 2. Experimental and simulated lu-29 scans across the GaN(0002) reflection 
of In^Gai-^N/GaN NW samples grown at different temperatures T a . The periodic 
In-cGai-^N QWs give rise to the first order superlattice peaks as indicated by the 
arrows. 



3. Structural analysis 

3.1. X-ray diffraction 

In order to determine the In concentration, we performed x-ray measurements with 
CuK«i radiation using a Ge(220) hybrid monochromator and a Ge(220) analyzer crystal. 
Symmetric Bragg scans across the GaN(0002) reflection are shown in figure |2j The 
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zeroth order satellite peak in the left shoulder of the GaN(0002) reflection indicates the 
increased average lattice spacing c avg of the MQW region due to In incorporation. This 
peak is surrounded by regularly spaced satellites on both sides, the ±l s ^ orders of which 
are marked by arrows in figure [2j These satellites arise from thickness interference at 
the periodic QW stack, and their positions can be used to determine the superlattice 
period d SL and c avg via Bragg's law [12] . We obtain <i SL = 19 ± 1 nm for these samples. 
For axial superlattices in NWs, lateral relaxation occurs. The XRD peak positions are 
determined by the average c-plane lattice parameters in the QW and barrier. For the 
interpretation of the XRD profiles, the following simplification may therefore be made: 
If the total height of the superlattice stack is larger than the NW diameter, the whole 
stack assumes a weighted average in-plane lattice constant [15]. In this case the Poisson 
effect vanishes for c avg [HI [15] . The average superlattice composition is hence given by 
Vegard's law: x avg = (c avg - c Ga N)/(ci n N - c Ga N)- 

3.2. TEM 

In the experiments shown in figure [2[ the contrast of the XRD satellite peaks is not 
sufficient to determine the QW thickness d we \\. From a calibration of the axial In^Gai-^N 
growth rate [12], we expect d we n = 5nm. To cross-check the structure of the NW 
samples, TEM was performed on the sample grown at T a = 608 °C. Figure [3] (a) shows 
a bright-field image of the sample cross section in two-beam condition with g = 0002 
at 200 kV. The Si substrate is located at the bottom. It can be seen that the NWs 
fluctuate in diameter, length and tilt. In the center of figure [3] (b) the complete 
heterostrucure of one individual NW is presented. The six-fold periodic structure of 
the axial In^Gai-^N/GaN MQW is visible. The single period highlighted by the white 
box is imaged by high-resolution TEM as depicted in figure [3] (c). Geometric phase 
analysis (GPA) reveals the In-induced strain in c-direction as shown in figure |3](d) [10J. 
The strain Ac/c is indicated with respect to the GaN barrier. It can be seen that it is 
reasonable to assume in the following the nominal value of d we ii = 5 nm. A line profile 
of Ac/c along the NW axis was obtained from the area in the GPA map highlighted 
by the black box, and is given in figure [3] (e). We would like to emphasize that due to 
the large ensemble fluctuations it is virtually impossible to obtain statistically relevant 
information about the average In content of the NW ensemble from such images. That 
is why our structural analysis focuses on techniques that probe the entire NW ensemble. 

3.3. X-ray simulation 

Using the method described in [Tl], we simulated the XRD profiles, and the resulting 
curves are also shown in figure El To account for the broadening of the satellite peaks, a 
standard deviation of the segment heights er^ = 2.5 nm has to be assumed. This rather 
large value reflects the inhomogeneity in heterostructure formation that is intrinsic to 
ensembles of self- induced nanowires grown by MBE [12]. For the present study, it is a 
decisive advantage that XRD probes a large number of NWs and directly reveals the 
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Figure 3. TEM analysis of axial In^Gai-^N QW in GaN NWs. (a) NW ensemble 
in bright field contrast, (b) NW tip with embedded MQW. (c) HRTEM image of one 
QW, from the area marked by the box in (b). (d) Ac/c map obtained from GPA, with 
the estimated QW thickness marked by dashed lines, (e) Ac/c strain profile along the 
NW axis, obtained within the black box in (d). The approximate location of the QW 
is marked red, and the QW strain estimate from XRD is given by the blue line. 



ensemble average and standard deviation. Compared to the QW thickness fluctuation 
within one sample of the series, the change in QW thickness between the different 
samples caused by different In concentrations is negligible. The composition of the 
QWs x is given by d we \\x = <isL^av g - Figure H] (a) shows x for d we n = 5.0 nm. x can be 
seen to decrease with increasing T a . The estimates for x resulting from the limiting cases 
for the QW thickness, d wc \\ ± ad, i.e. between 2.5 nm and 7.5 nm, are shown as error 
bars. For the sample grown at T a = 608 °C the strain Ac/c, as obtained from the XRD 
data, is marked for comparison in the strain profile from GPA in figure |3] (e) . The solid 
blue line corresponds to d wc n = 5.0nm, Ac/c = 1.75% and x = 11.5%. The dashed 
lines give the estimates from XRD for the limiting QW thickness estimates as discussed 
above. Within the large experimental error, the data agree between this individual NW 
in HRTEM and the NW ensemble in XRD. 
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Figure 4. In content x in In^Gai-^N/GaN NW axial superlattices grown at different 
temperatures T a . (a) x from XRD satellite peak positions, (b) x derived from E'i(LO) 
phonon frequency in the Raman spectra. The uncertainty ranges are described in the 
main text. 




Figure 5. Raman spectra excited at 3.0 eV for In^Gai-^N/GaN nanowire axial 
superlattices grown at different temperatures T a . The dashed line marks the i?i(LO) 
phonon frequency peak position of the sample with the highest In content. 



3.4- Raman spectroscopy 

A complementary estimation of x was achieved by resonant Raman spectroscopy [16J. 
The frequency of LO phonons in In^Gai-^N depends on both the In content and strain 
|17j . In order to selectively enhance the Raman signal from the In x Gai_ x N QWs, we 
excited the spectra with a photon energy of 3.0 eV, i. e., close to the resonance with the 
fundamental band gap for x ~ 10 %. The efficiency of scattering by LO phonons via 
the Frohlich mechanism can be considerably enhanced by choosing this kind of resonant 
Raman spectroscopy [18]. In fact, the resonantly enhanced LO-phonon signal from the 
In^Gai-^N QWs is much stronger than that from the GaN base material. Using non- 
resonant conditions, Raman spectra from In^Gai-^N would not be detectable because of 
the small scattering volume provided by the thin QWs and the small separation between 
In^Gai-^N and GaN phonon frequencies at low In contents. The corresponding spectra 
of our samples are depicted in figure [5) They are dominated by the -Ei(LO) phonon 
peak. The shoulder at 685 cm -1 might be attributed to the S band which is associated 
in literature with disorder or defects [TTJ EHj. Contrary to planar layers, light is coupled 
into and out of c-plane In^Gai-^N laterally through the side facets (perp. to the c- 
axis) in our as-grown ensembles of vertical NWs, and scattering by E'i(LO) phonons 
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Figure 6. Room temperature PL spectra of In^Gai-^N/GaN NW axial supcrlattices 
grown at different temperatures T a . (a) Excitation at 3.8 eV (logarithmic), (b) Direct 
excitation at 3.0 eV (linear). 

dominates over that by Ai(LO) phonons [2D] . Since the -Ei(LO) frequency shift with In 
content and strain is not well documented, we assumed that the offset against yL(LO) 
phonons is independent of the In content and strain, and used published dependencies 
for A\{LO). With decreasing T a , the peak shifts towards lower frequencies. This shift 
amounts to — 149xcm~ 1 for relaxed In x Gai_ x N [ITl [21] . Biaxial compressive strain 
induces higher LO phonon frequencies, and we derive a shift of — 58a; cm -1 for In x Gai_ x N 
grown pseudomorphically on GaN from the data given in [22] for GaN. These limiting 
cases are indicated by the error bars in figure H] (b). The values of x plotted in this 
diagram were derived under the assumption that the active region relaxes laterally and 
the relative strain is equal to (ibarricr/^SL as extracted from the XRD measurements. In 
agreement with the XRD results, x decreases with increasing T a . 

The range of In content deduced from the Raman measurements is smaller than the 
range derived from XRD. This difference could arise from the resonant excitation which 
favors regions with x ~ 10 % if the composition fluctuates within a given sample. The 
weaker Raman signal for samples with higher T a as seen in figure [5] is then explained by 
the detuning of the resonance condition for lower In content. Another factor might be the 
superimposed contribution of ^(LO) phonons, which would lead to an overestimation of 
x. While for the rest of this discussion we rely on the absolute values of x as derived from 
XRD, we do see the resonant Raman results as an independent proof of the systematic 
In content variation with growth temperature. 

4. Luminescence properties 

4-1. Photoluminescence 

With the structural parameters elucidated, we turn to the emission properties. Room 
temperature photoluminescence (PL) spectra of the different NW ensembles are 
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Figure 7. Experimental PL peak positions and calculated planar In^Gai-^N QW 
transition energies. The shaded area results from the maximum and minimum 
estimation of x and dwell from XRD. 



presented in figure [6] (a) for excitation at 3.8 eV, i.e. above the band gap of GaN. 
They show three main emission bands. The peak at 3.4 eV represents the free exciton 
transition in the strain-free GaN NW base and is common to all samples. The spectra 
are normalized to this peak. The second peak at 3.0 eV, which is observed for all samples 
grown at lower temperatures, is attributed to emission from insertions of cubic GaN (c- 
GaN) [23]. Direct evidence for the formation of c-GaN was obtained in situ by RHEED 
as described above. The redshift with respect to the c-GaN room temperature band 
edge of 3.2 eV may be explained by the quantum-confined Stark effect (QCSE) in thin 
layers of c-GaN in a matrix of hexagonal GaN [24] . The third PL emission peak in the 
energy range from 2.2 eV to 2.5 eV shifts systematically with T a . In view of the variation 
in In content with T a , we attribute this emission to the In^Gai^N QWs. The shift in 
the peak position can be seen more clearly in the normalized spectra depicted in figure [6] 
(b) which were obtained by direct excitation at 3.0 eV. The intensity of the In^Gai-^N 
emission decreases drastically with decreasing In content. This decrease in intensity with 
decreasing In content (and therefore presumably higher crystalline quality and reduced 
QCSE) is unexpected and contradictory to experience with planar (In,Ga)N structures. 

4-2. Calculation of transition energies 

Figure [7J shows the observed peak positions for the emission from c-GaN and the 
In^Gai-^N QWs as a function of T a . To correlate the experimentally observed emission 
energies and structural parameters, we calculated the expected transition energies by 
solving the one-dimensional Poisson and Schrodinger equations self-consistently [25], 
As input, we used the In^Gai-^N QW thicknesses G? we ii given above, compositions x 
derived from XRD, and published values for the bowed band-gaps [26], the band- 
offset [27], and the bowed polarization [28J. The simulated transition energies are 
added to figure [7J for comparison with the experimentally observed PL peak positions. 
The shaded range of values results from the estimates for x derived from XRD with 
different assumptions for <i we n as shown in figure H] (a). The PL at 3.0 eV is of almost 
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constant energy, which supports our assignment of this emission to c-GaN insertions. 
In contrast, the luminescence around 2.2-2.5 eV, assigned to emission from In r Gai_. ! ;N, 
systematically shifts to higher energies with decreasing In content. However, this shift is 
significantly smaller than the calculation predicts, regardless of the specific assumption 
for dweii. Thus, the situation in our NW based axial In x Gai_ x N QWs differs considerably 
from that of planar QWs, for which one- dimensional calculations are in much better 
agreement with the experimental emission energies [29]. We have recently shown that 
the reduction in piezoelectric polarization from better strain accommodation in the NW 
structure is not significant enough to explain the reduced shift in emission energy [30] . 

5. Discussion 

Both the luminescence energy range we observe and the absence of bright emission 
above 2.7 eV normally expected for low In content are in accord with similar optical 
measurements in other reports where the In content was not determined independently 
[H El El HI El El Ej. This agreement is intriguing because it indicates a universal 
characteristic of axial In^Gai-^N/GaN NW heterostructures. In principle, low PL 
intensity could be caused by nonradiative point defects, but it seems highly unlikely 
that the density of these hypothetical point defects increases with decreasing In content 
at higher growth temperatures. Moreover, point defects would not affect the emission 
energy of the In x Gai_ x N/GaN heterostructures. As a more probable explanation, three 
types of carrier localization have been discussed to influence the emission energy and 
intensity of such heterostructures. First, exciton localization induced by compositional 
fluctuations was observed in similar samples [91 [301 EH]- Spatially direct exciton 
localization alone, however, would explain neither the low absolute intensity nor the 
cross-over between the observed PL energy and calculated transition energy for the 
highest In content seen in figure [71 Second, we have recently seen that individual 
electron and hole localization in separate potential minima with varying distance plays 
a role for In x .Ga 1 _ x N QWs in GaN NWs [30]. A spatial separation of electrons and 
holes implies a reduced oscillator strength, thus accounting for the low PL intensity. 
If this separation occurs mainly in the lateral direction which is not affected by the 
QCSE, the redshift could be reduced compared with the one-dimensional calculations 
in axial direction. Third, a three-dimensional numerical analysis [32] indicates that axial 
In^Gai-^N/GaN NW heterostructures exhibit laterally inhomogeneous strain resulting 
in band bending and a lateral separation of electrons and holes. Likewise, this could 
explain the low PL intensity. Moreover, numerical studies for axial GaN/(Al,Ga)N 
QWs in NWs demonstrate that both geometry and chemical composition significantly 
affect the location of the potential minima in the conduction and valence band leading 
to carrier localization [33J El]- Further, our TEM images as well as those found in 
the literature [U El [10] show that the In x Gai_ x N QWs are laterally embedded in a 
GaN shell in contrast to the structure assumed in [32]. Thus, although the detailed 
findings of [32J on where electrons and holes are localized cannot be generalized, the 
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inhomogeneous strain could possibly account for the cross-over between the results of 
the one-dimensional model and the observed PL. 

In conclusion, our results show a clear trend of emission energy with In 
incorporation, even though it is weaker than expected. In contrast to planar 
In x Gai_ x N/GaN heterostructures, the analysis of the optical transitions in NW 
heterostructures necessitates — even for a rather basic level — fully three-dimensional 
simulations that take into account lateral variations in strain. Our findings 
suggest, however, that the interplay of several mechanisms for carrier localization 
has to be considered to describe the luminescence from NW based Inj.Gai_ x N/GaN 
heterostructures. 
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